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ABSTRACT 

Polycyclic aromatic hydrocarbons (PAHs) have been proposed as potential 
carriers of the unidentified infrared bands (UIRs) and the diffuse interstellar 
bands (DIBs). PAHs are not likely to form by gas-phase or solid-state interstellar 
chemistry, but rather might be produced in the outflows of carbon-rich evolved 
stars. PAHs could form from acetylene addition to the phenyl radical (CgHs), 
which is closely chemically related to benzene (CeHg) and ort/io-benzyne (o- 
C6H4) . To date, circumstellar chemical models have been limited to only a partial 
treatment of benzene-related chemistry, and so the expected abundances of these 
species are unclear. A detection of benzene has been reported in the envelope 
of the proto-planetary nebula (PPN) CRL 618, but no other benzene-related 
species has been detected in this or any other source. The spectrum of 0-C6H4 is 
signiflcantly simpler and stronger than that of CeHs, and so we conducted deep 
Ku-, K- and Q-band searches for 0-C6H4 with the Green Bank Telescope. No 
transitions were detected, but an upper limit on the column density of 8.4x10^^ 
cm~^ has been determined. This limit can be used to constrain chemical models 
of PPNe, and this study illustrates the need for complete revision of these models 
to include the full set of benzene-related chemistry. 
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Subject headings: astrochemistry circumstellar matter - stars: individual (CRL 
618) radio lines: stars 



Introduction 



Polycyclic aromatic hydrocarbons (PAHs) are large molecules with carbon atoms ar- 
ranged in five- or six-membered rings and are thought to form from smaller carbon ring 
molecules such as benzene (CeHg) and its derivatives. PAHs are very stable against pho- 
todissociation, and so could be present in interstellar and/or circumstellar environments. The 
unidentified infrared bands (UIRs) occur at frequencies characteristic of aromatic molecules, 
and so PAHs have been suggested as potential UIR carriers. Likewise, PAHs have been 
proposed as carriers of the diff use int e rstell ar bands (DIBs) . We refer the reader to a recent 
review of the DIB problem by ISarre I (120061 ) and references therein for a discussion of PAHs 
and their relation to the UIRs and DIBs. Significant laboratory spectroscopic work has been 
dedicated to PAHs and other benzene-related species to support astronomical observations, 
yet none of these species has been unambiguously detected in space. Only benzene has been 



tentat ively detected toward the proto-planetary nebula (PPN) CRL 618 (jCernicharo et al 



2001al ). One vibrational band of benzene was observed, and a column density of 5x10 



15 



cm ^ and a kinetic temperature of 200 K were determined. 



CRL 618 is an evolved PPN with a centr al BO star and an ultra-compact H II region sur 



rounded by a carbon-rich molecular envelope (jCernicharo et al. 



20041 ). This source has optical high velocity bipolar out flows (jTrammell 1 120001) in additio n 



to a low velocity expanding torus of molecular emission (j Sanchez Contreras fc Sahai 



2001b; Sanchez Contreras fc Sahai 



CRL 6 18 has a rich molecular inventory including a variety of hydrocarbons [see lPardo et al. 



2004). 



(120071 ) and references th erein], but searches fo r biologically important molecules have pro- 
vided only upper limits (IRemijan et al. 1120051 ). 



Very few observational constraints have been placed on the chemical and physical prop- 
erties of PPNe, especially regarding benzene-related chemistry. Observational studies of the 
simplest benzene derivatives are the essential flrst steps to understanding the formation of 
PAHs in PPNe. The carbon-rich nature of CRL 618, coupled with the benzene detection, 
makes it an excellent target for initial benzene derivative searches. We have therefore con- 
ducted a search for the benzene derivative ortho-henzjne (0-C6H4) in CRL 618 with the 
Green Bank Telescope (GBT) Ku-, K-, and Q-band receivers. An overview of circumstel- 
lar benzene-related chemistry, details of the observations, the results of our search, and a 
discussion of these results are presented below. 
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2. Benzene-Related Circumstellar Chemistry 



PAHs have long been thought to be produced in carbon-star outflows, as it seems 
unlikely that such large molecules could be produced by gas-phase or grain-assisted chemistry 
within the interstellar medium. Although the PPN phase of stellar evolution is typically 
very short, lasting only about ten thousand years, the chemistry that takes place in the 
envelope of the post-asymptotic giant branch (AGB) star leads to the formation of complex 
hydrocarbons. PPN chemistry is steeply dependent on density because circumstellar material 
undergoes energetic processing, and formation mechanisms for large molecules are likely 
a combination of radical-molecule and ion-molecule reactions. This is in contrast to low 
density interstellar chemistry that is almost exclusively driven by ion-molecule reactions. 
Both radical-molecule and ion-molecule reactions can lead to the formation of benzene and 
related species, and su c h benzene formation mechanisms have been included in PPN models 



( iRedman et al. I l2003l : IWoods et al. I l2003l ). However, no PPN model includes both the 



radical- and ion-based benzene reaction networks. This striking deflciency in PPN models 
makes interpretation of observational results quite difficult, and the radical- and ion-based 
chemical mechanisms must be considered separately until more complete PPN models are 
developed. 

The torus of a PPN is a high energy, high density en vironment, and the chemistr y 



is expected to be similar to that observed in combustion (IFrenklach fc Feigelson I Il989l ). 



Radical-based mechanisms based on combustion chemistry have been proposed for b enzene 



and PAH formation in the circumstellar shells of AGB stars (jCherchneff et al. Ill992l ). and a 



subset of these reactio ns was also used to model the chemistry of clumps during PN evolution 
( IRedman et al. II2003I ). A summary of radical-driven benzene chemistry is shown in Figure 



[U Here, hydrocarbon radicals react with acetylene (C2H2), and ring-closure forms 0-C6H4. 
Additional reactions lead to the phenyl radical (CeHs) and benzene. Other routes to benzene 
include either reaction of two C3H3 radicals or, alternatively, the formation of straight-chain 
energetic CqRy (n-CgHy), which can then undergo ring-closure. 

As a PPN evolves and its surrounding gas expands, the density decreases and the 
material in the circumstellar envelope is subjected to photoprocessing from the central star, 
leading to ion-molec ule chemistry. An ion -molecule benzene formation mechanism has been 
proposed for PPNe (IWoods et al. II2003I ). and a summary of this mechanism is shown in 
Figure [21 In this network, hydrocarbon ions react with acetylene, and ring closure forms 
CeH^. Additional hydrogenation and/or electron recombination leads to 0-C6H4, c-CeH^" , 
and ultimately CeHg. It should be noted that the original PPN model investigated only 
the prim ary route to benzene shown in Figure [2] and did not include the 0-C6H4 formation 



reaction (IWoods et al. II2003I ). 
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One of the benzene derivatives formed in these reactions, CfiH^, is a suggested precur- 



sor to PAHs in circumstellar environments (IFrenklach fc Feigelson I Il989l : ICherchneff et al 



19921 ). As is illustrated in Figure [31 CgHs can undergo a series of radical- molecule reactions 
involving acetylene to produce a naphthalene-like species (CioHg) via ring closure. Subse- 
quent reactions of this nature can lead to larger PAHs. 

Studies of other high energy environments indicate that CgHs and 0-C6H4 may well 
coexist with benzene in CRL 618 if circu mstellar chemistry is si r nilar to combustion or 
plasma chemistry. The reaction network of IFrenklach fc Feigelson I (119891 ) is based on soot 
production mechanisms in hydrocarbon flames, and theoretical studies show that benzene 



unim olecular decomposition leads to CeHs and 0-C6H4 during co mbustion (|Mebel et al 



2OOII ). Benzene electrical discharges also produce CqB.^ and 0-C6H4 (IMcMahon et al. II2003I ). 
There are no reliable predictions of CqHq, CqB.^, and 0-C6H4 abundances in PPNe because 
of the partial treatment of their chemistry in models, but observations of CeHs and 0-C6H4 
would provide important limits for future modelling studies. We therefore began a search 
for these species in CRL 618. 

CgHg has no permanent dipole moment and can therefore only be studied in the in- 
frared, but both CeHs and 0-C6H4 can be probed by radioastronomical tec hniques. Th e 



rotational spectra of o-( 


-1 TT 

^fiil4 


1986 




Robertson et al. 


2003 



20031 : iKukohch et al. II2003I : iMcMahon et al. 1120031 ). Both species 



are carbon ring, planar, near-oblate asymmetric rotors with C2,, symmetry along their b 
inertial axes. The calculat ed dipole moment of o -Cp,B.4 is 1.38 D (IKraka fc Cremer Ill993l ) 
and that of CeHs is 0.9 D (IMcMahon et al. 1120031). The unpaired electro n in CqB.^ leads to 
hyperfine splitting of the lower rotational states (IMcMahon et al. II2003I ). and so there are 
many rotational transitions to sample observationally. However, the larger dipole moment 
coupled with the lack of hyperfine splitting yields much stronger lines for 0-C6H4 and makes 
this species a more likely candidate for detection if the abundances are similar. 



3. Observations 

Observations of 0-C6H4 were conducted with the NRAoEl 100 m Robert C. Byrd GBT 
between 2006 September 4 - 2007 January 28 using the Ku-band (12 - 15.4 GHz), K-band 
(18 - 22.5 and 22 - 26.5 GHz,) and Q-band (40 - 48 GHz) receivers. The eight intermediate- 
frequency (IF), 200 MHz, three- level GBT spectrometer configuration mode was used, pro- 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation, operated 
under cooperative agreement by Associated Universities, Inc. 
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viding four 200 MHz frequency bands with a channel separation of 24.4 kHz in two polar- 
izations. The assumed CRL 618 J2000 pointing position and LSR source velocity were a= 
04''42™'53*.7, 5 = -|-36°06'53".0 and -27.5 km s~^, respectively. Data were acquired in the 
OFF-ON position-switching mode. A scan included two-minute integrations for each posi- 
tion beginning with the OFF-source position, which was located 60' East in azimuth of the 
ON-source positi on. Antenna temper atures with estimated 20% uncertainties were recorded 
on the T\ scale (jUlich fc Haas Ill976l ). The GET half-power beamwidths are given by Q\, = 
740"/z/(GHz). Dynamic pointing and focusing corrections were applied and observations of 
the quasar 0359+509 were used to adjust the zero points every two hours or less. Observa- 
tions from multiple nights and both polarizations were averaged for each frequency window, 
and the data were Hanning smoothed over three channels with the GBDish data reduction 
program. The resultant Q-band spectrum is shown in Figure HI 

Table 1 lists the 0-C6H4 rotational transitions in the observed frequency windows. The 
transition quantum numbers, transition rest frequency (z/), the Einstein A coefficient of the 
transition times the upper state degeneracy (Agu), the energy of the upper level (Eu), the 
observed l-cx RMS level [Tmb), and the beam efficiency (77) are listed in the first six columns. 
The transitio n frequencies and in tensities are from the Cologne Database for Molecular 
Spectroscopy (iMiiller et al. II2005I ). The d.rms routine in GBDish was used to calculate the 
RMS level in the scale for each spectral window. The RMS level is equal to the standard 
deviation of the noise in a line-free region and was calculated after the Hanning smoothing 
was applied. These values were then converted to t he Tmr scale by the relationship Tmb = 
T^/rj. The value o f r) was derived from a fit to the iRuze I (119661 ) formulation, as is outlined 
in Equation (2) of IHoUIs et a"l~l J2007h . 



4. Results and Discussion 

No spectral features associated with 0-C6H4 were detected during this search, and there 
are no unidentified spectral features in any passband. Figure S] shows a plot of the observed 
CRL 618 Q-band spectrum overlaid by a predicted 0-C6H4 spectrum at a column density of 
10^^ cm~^, a rotational temperature of 200 K, and a source size of 10". If these parameters 
are representative of 0-C6H4, the blended 60,6 — 5i,5 and 61,6 — 5o,5 transitions at 40828.1686 
and 40829.9929 MHz, respectively, should have been easily detected at the 30 mK level. 

The RMS levels reached during the observations allow calculation of the 0-C6H4 column 
density upper limit in CRL 618, and the 3-a upper limits determined from each observed 
transition are presented in Table [H These column density upper limits were calculated using 
the following expression, adapted from equation (1) of Nummelin et al. (1998): 
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where Nt is the beam averaged total column density, T},dv is the transition integrated 
intensity, Q{Trot) is the rotational partition function, and T^ot is the molecular rotational 
temperature. Since no lines were observed, the integrated intensity was approximated as 
that of a line with a peak intensity at the 3-0" RMS level, Trmsi aiid an assumed linewidth 
Af. The value of Trms was determined by Trms=^Tmb/ where n is the number of 
channels across the linewidth Af. The factor does not include a complete statistical 

treatment of the assumed Gaussian lineshape, but it does approximate the overestimation of 
the RMS level from spectral oversampling. Examination of previously reported upper limit 
calculations indicates that, in nearly all cases, this factor is either neglected entirely, or at 
the very least not discussed explicitly. We find this to be a gross omission for instruments 
such as the GBT spectrometer, where the channel width is often significantly smaller than 
the linewidth. From the Trms values, correction for beam dilution gave Tf, through the 
expression Tb=TjiMs/B, where the beam filling factor, B, was calculated from the source 
size, 9s, and the beam size, 9b, by the relationship B = 9'^/[9'^ + 9l]. The integrated intensity 
was then approximated as J^^Tb dv = 1.064 T;, Av, where the 1.064 factor arises from the 
assumed Gaussian lineshape. We have included this factor for completeness, although it is 
likely minimal compared to the uncertainty in T^. 

The calculated 3-a 0-C6H4 upper limits for CRL 618 are presented in Table [H These 
calculations required assumptions for the values of Trot, , and 9s. The linewidth was 
assumed to be that observed for other identified lines, 10 km s~^. Determination of an 
appropriate temperature and source size was less straightforward, as there is much discrep- 
ancy in t he literature regarding th ese values. The kinetic temperature derived for benzene 
is 200 K (jCernicharo et al~ 2001a ). but an IRAM 30- meter millimeter line su rvey indicates 



temperatures of 250 - 275 K for the torus and 30 K for the circumstellar shell (jPardo et al. 



20071 ). It is expected that 0-C6H4 would have very similar physical properties to those of 
benzene in this source, and so a temperature of 200 K was assumed for the upper limit 
calculations. 



There is similarly contradicting information regarding the CRL 618 source size. iPardo et al. 



(120071 ) found torus and circumstellar shell sizes of 1.5" and 3.0 - 4.5", respectively, while mil- 



limeter interferometric observations with a ~5" beam gaye sou rce sizes <10" for molecules 



expected to trace the extended envelope (IRemijan et al. II2005I ). A molecule such as 0-C6H4 



is likely to be present in the torus where higher densities shield it from photodissociation, 
but such a species could also be present in more extended regions if readily formed by ion- 
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molecule chemistry. The GBT beam is >16", and so our observations probed both the torus 
and the extended envelope of CRL 618. A source size of 10" was assumed for the upper hmit 
calculations, as this includes the entire molecular envelope and therefore all possible regions 
of emission. 

As is shown in Table [H the most strict 0-C6H4 column density upper limit determined 
from t hese observations is 8.4x 10^^ c m~^. The observed benz ene column density is 5 x 10^^ 



cm ^ (jCernicharo et al. Il2001al ). and lMcMahon et al. I (120031 ) report a CgHs column density 



upper limit in CRL 618 of 4 xlO^^ cm~^. The 0-C6H4 upper limit is therefore much lower 
than the limits for the related species. Given the incomplete nature of the chemical models, 
however, the significance of this limit is unclear. 

Further interpretation of this limit will require additional observational and modelling 
studies, as this work highlights the nearly total lack of information regarding the physical 
and chemical properties relevant to benzene-related chemistry in circumstellar environments. 
Interferometric observations would prove quite useful, as such studies would probe the spatial 
distribution of molecules in the torus and extended envelope, eliminating the uncertainty 
in beam dilution effects. Laboratory measurements of the millimeter and submillimeter 
spectra of 0-C6H4 would aid these observations. We also strongly encourage complete revision 
of existing PPN models to investigate both the radical- and ion-based benzene chemical 
networks for comparison to 0-C6H4, CgHs, and CqB.q observations. Only a combination of 
further observations, modelling, and laboratory measurements of benzene derivatives will 
lead to full understanding of PAH formation mechanisms in PPNe. 
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Table 1: Summary of 0-C6H4 Observations Toward CRL 618 



Nt 

'^'k'^,k'^~'^k^,k'J ^ Eu Tmb V Upper Limit 

(MHz) (s-^) (K) (mK) (cm-^) 



4 


2,2 


-4i,3 


13250.8230 


6.08x10 


3 


3,1 


— 32,2 


13680.0050 


3.59x10 


6 


5,1 


- 64,2 


14608.9021 


1.28x10 


3 


1,2 


— 3o,3 


14619.6969 


2.40x10 


2 


0,2 


-ll,l 


14753.9974 


4.23x10 


3 


0,3 


-2l,2 


21770.7375 


1.60x10 


4 


2,3 


-4i,4 


21963.3591 


7.93x10 


11 


8,3 ■ 


- 117,4 


22142.0552 


5.35x10 


9 


5,4 


-94,5 


22143.1890 


4.10x10 


3 


1,3 


— 2o,2 


22216.3288 


2.88x10 


7 


5,3 


— 74,4 


23578.9636 


5.35x10 


2 


2,1 


-li,o 


24109.5097 


1.40x10 


3 


1,2 


-22,1 


24842.4030 


8.22x10 


8 


2,6 


— 81,7 


40403.3815 


1.64x10 


8 


3,6 


— 82,7 


40421.4144 


9.83x10 


5 


1,4 


-42,3 


40595.9260 


1.40x10 


22i3,9 - 


- 22i2,io 


40649.5538 


1.13x10 


7 


1,6 


— 7o,7 


40674.4124 


5.21x10 


7 


2,6 


-7i,7 


40677.1133 


8.69x10 


6 


0,6 


- 5l,5 


40828.1686 


4.05x10 


6 


1,6 


- 5o,5 


40829.9929 


2.43x10 


3 


3,0 


-22,1 


41071.9441 


8.57x10 



■7 


5.28 


1.48 


0.92 


3.47x10^^ 


■7 


3.70 


1.76 


0.92 


6.87x10^^ 


6 


12.37 


1.45 


0.91 


1.59x10^^ 


■7 


2.96 


1.45 


0.91 


8.22x10^5 


7 


1.19 


1.45 


0.91 


4.62x10^5 


6 


2.26 


2.47 


0.86 


1.79x10^^ 


7 


4.67 


2.47 


0.86 


3.63x10^^ 


6 


37.44 


2.18 


0.85 


5.28x10^4 


6 


24.25 


2.18 


0.85 


6.60x10^^ 


6 


2.26 


2.18 


0.85 


8.72x10^'' 


6 


15.55 


2.51 


0.84 


5.56x10^'' 


6 


1.77 


1.76 


0.84 


1.42x10^^ 


■7 


2.96 


2.13 


0.83 


2.90x10^^ 


5 


16.25 


6.42 


0.64 


4.19x10^4 


6 


16.25 


6.42 


0.64 


6.98x10^^ 


■5 


6.62 


5.92 


0.63 


4.37x10^"^ 


■4 


137.89 


5.92 


0.63 


8.40x10^3 


6 


11.45 


5.92 


0.63 


1.19x10^^ 


6 


11.45 


5.92 


0.63 


7.16x10^^ 


5 


7.23 


6.21 


0.63 


1.59x10^^ 


■5 


7.23 


6.21 


0.63 


2.64x10^^ 


6 


3.74 


9.40 


0.63 


1.12x10^^ 
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Fig. 1. — Radical-based reaction scheme for PPNe benzene chemistry based on 
Frenklach &: Feigelson I (119891 ). 
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Fig. 2. — Ion-based chemical scheme for PPNe benzene chemistry based on IWoods et ah 
J2003h . 
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Fig. 3. — Potential PAH for mation routes in circumstellar shells (IFrenklach &: Feigelson 
19891 : ICherchneff et al. lEoQ^ ). 
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Fig. 4. — Observed CRL 618 Q-band spectrum (grey) overlaid with a simulation of the 
0-C6H4 spectrum (black) at a column density of 10^^ cm~^, a rotational temperature of 200 
K, and a source size of 10". The observed spectrum is comprised of four 200 MHz windows 
and has been smoothed to a resolution of 200 kHz. The features at 40500 and 40900 MHz 
arise from noise at the edges of the spectral windows. 



